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Abstract 



We demonstrate that, for every grand unified model based on a conven- 
tional type I seesaw mechanism leading to a normal light neutrino mass hi- 
erarchy, one can easily generate a corresponding model with an inverted hi- 
erarchy which yields the same neutrino oscillation parameters. However, the 
latter type model has several unattractive instabilities which will disfavor any 
grand unified type I seesaw model, if an inverted neutrino mass hierarchy is 
observed experimentally. This should be contrasted with the softly-broken 
Lg — — L-T fiavor symmetry models which are eliminated, if the data favors 
a normal mass hierarchy. 
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Since the first convincing evidence in 1998 for atmospheric neutrino oscillations by the 
Super-Kamiokande (S-K) collaboration and subsequent confirmation of hints of solar neu- 
trino oscillations by both the Sudbury Neutrino Observatory and KamLAND collaborations, 
great progress has been made in pinning down the values of the leading oscillation parame- 
ters Am32, sin^ 2^23, Am|i and sin^ 2^12 [1]. The same can not be said for the other 3-family 
mixing parameters, sin^ 26i^, and the leptonic Dirac and Majorana CP phases 5, Xi ^"^^ Xi- 
The angle ^^13 has been bounded by the CHOOZ and Palo Verde reactor experiments [2] 
to be noticeably smaller than the other two angles, while the three phases are completely 
unknown. Moreover, most important for the subject of this letter, even the 3-neutrino mass 
hierarchy, normal or inverted, remains undetermined. Experimental determination of the 
hierarchy and ^13 will be critical for the continued viability of models in the literature which 
attempt to explain the neutrino mixing data; cf. [3] for extensive reviews of such models. 

For models based on some flavor symmetry, those exhibiting a softly broken L^ — L^ — L^. 
lepton number produce only an inverted hierarchy [4]. In contrast, all of the successful grand 
unified models in the literature employing a conventional type I seesaw mechanism appear 
to prefer a normal hierarchy; cf. [5] for some still viable models of this type. We seek to 
explore here why this is the case. With such a dichotomy present, clearly identification of 
the correct hierarchy observed in Nature will serve to eliminate a large batch of models of 
one kind or the other. On the other hand, grand unified models based on a type II [6] or 
type III [7] seesaw with direct or induced Higgs triplet contributions can be constructed by 
design to have either a normal or an inverted hierarchy. 

In this letter we show how, given a satisfactory grand unified model with normal neutrino 
mass hierarchy based on a type I seesaw, one can easily generate a corresponding model with 
an inverted hierarchy which also satisfies the same known neutrino oscillation parameters. 
The issue becomes which hierarchy model is more satisfactory theoretically in terms of flavor 
symmetry, heavy neutrino mass hierarchy with possibly successful leptogenesis, and stability 
of the solution. To address this point, we use a much studied 5*0(10) grand unified model 
with a ?7(1) X Z2 X Z2 flavor symmetry proposed by the author in collaboration with S.M. 
Barr which leads naturally to a normal hierarchy [8]. Results for the light neutrino mass 
matrices are also derived for the class of models in which the charged lepton mass matrix 
is diagonal in flavor space. Features obtained in the models considered are more generally 
expected to hold for the whole class of grand unified models with a type I seesaw mechanism. 

We begin with the known neutrino oscillation data and the definitions that apply for the 
parameters appropriate for a normal vs. inverted hierarchy. Experimentally the atmospheric 
and solar neutrino mixing parameters are approximately equal to [1,2] 

|Am|2| ^ 2.5 X 10-3 eV^, sin2 2^23 ^ 1-0, 

Amii ~ 7.0 X 10-5 ev2, tan^ 612 ^ 0.40, (1) 

Am|i ~ Ami2, sm^2ei^ < 0.16. 

Although the mass squared differences are determined, the actual neutrino mass scale is 
undetermined as is the mass hierarchy, i.e., whether Am^j = m\ — m\ is positive or negative. 
Clearly, however, the larger (smaller) mass separation is associated with the atmospheric 
(solar) neutrino fiavor mixing. For the solar neutrinos, the mass squared difference Am2i is 
known to be positive, since the solar mixing angle ^12 lies in the first octant. 
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The Maki-Nakagawa-Sakata (MNS) neutrino mixing matrix with elements Uai relates 
the flavor state a = e, ji, r to the mass eigenstates i = 1, 2, 3 in their respective orders. 
By convention, we have 



/ C12C13 S12C13 sise 

Umns - 



— S12C23 — Ci2S23Si3e*'' C12C23 — Si2S23Si3e*'^ S23C13 
V S23S12 - Cl2C23Sl3e*'^ -C12S23 - Si2C23Sl3e*'' C23C13 / 



(2) 

where the cosines and sines refer to rotations of the mass eigenstates. From the mixing data 
in Eq. (1), we can approximate C23 — S23 — 1/a/2, su — ^2/7, cu — ^5/7, and S13 ^ 0.2. 
The mass eigenstates are then approximately given in terms of the flavor eigenstates by 

i/3 ~ Si3e"'i/e + \/ici3(i/^ + i/^), 

where the states v\ and V2 are involved with the solar neutrino transitions, while V2 and 2/3 are 
involved with the atmospheric neutrino transitions. As such, the same Umns mixing matrix 
leading to the same mass eigenstate flavor compositions applies, whether the hierarchy is 
normal, < 1712 < m^, or inverted, m'g < m'^ < m'2, as illustrated in Fig. 1. For the 
inverted case, Am^^ = — m'^ and |Am|2| = Tn'2 — with the mass square differences 
specified experimentally in Eq. (1). 

We now turn to the issue of inverted vs. normal mass hierarchy in grand unified models 
with the conventional type I seesaw mechanism. Since most, if not all, such models in the 
literature [5] favor a normal hierarchy, we start with that case. The specified Dirac and 
right-handed Majorana mass matrices, Mjv and Mr, will both have normal hierarchy in 
that their mass eigenvalues with the smaller values will be more closely spaced. The light 
left-handed neutrino mass matrix follows from the conventional type I seesaw formula [9] 
according to 

= -MnM^'mI, (4) 

which by assumption also has a normal hierarchy, albeit a rather mild one. This mass matrix 
is diagonalized by the unitary transformation, U^^, whereby 

Ul^MlM,U,^^^i&g{ml, ml, ml). (5) 

Since the Mj, matrix is complex symmetric, one can choose phases for U^j^ such that diago- 
nalization of M,^ itself occurs with 

C/jM^C/^^ = diag(mi, m2, m3), (6) 

where the mass eigenvalues are real and positive with mi < m2 < m^. The MNS mixing 
matrix is given by the product of the two unitary matrices, Ul^ which diagonalizes the 
charged lepton mass matrix, and U,^^ above according to 
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Vmns = Uj^JJi^j^ — Umns^i 
$ = diag(e^>^i, e^^^, 1), 



(7) 



where Umns is specified by the convention of Eq. (2), and $ is a diagonal phase matrix 
involving the two Majorana phases, xi ^"^^ X2- 

Note that if the right-handed Majorana mass matrix had the structure Mr oc MJ^Mn, 
by the seesaw formula in Eq. (4), M^^ would be proportional to the identity matrix with 
three-fold degeneracy. This suggests that some satisfactory may be found leading to 
an inverted neutrino mass hierarchy. We now determine the corresponding matrices which 
yield not only the correct mass spacing with an inverted hierarchy, but also the correct MNS 
mixing matrix. 

For this purpose, let us choose the heaviest light neutrino mass m.2 to be assigned some 
value to be specified later. The other light neutrino masses can then be determined from 
the observed mass squared differences in Eq. (1) and m'2 according to 



We invert the counterpart of Eq. (6) written for the inverted hierarchy case and identify 
the resultant M'^ with the seesaw formula involving the same Dirac neutrino matrix but the 
new right-handed Majorana matrix, Mj^: 



By using the same unitary tranformation 17,^^ for the inverted hierarchy case as in the normal 
hierarchy case and the masses with 777,3 < '^i < ^2 ordered in the diagonal matrix as 
above, we guarantee that the same desired MNS mixing matrix, Umns, a-nd Majorana phase 
matrix, $, are obtained in the inverted case. In general one will find that the mild normal 
hierarchy for the light neutrino mass matrix, M^,, is replaced by a completely different 
texture with no simple flavor symmetry obvious for the inverted hierarchy version. The new 
right-handed Majorana matrix which accomplishes this follows from the seesaw formula in 
Eq. (9) by inversion: 



The procedure described above for determining the new light neutrino mass matrix 
which obtains in the inverted hierarchy case, given the original normal hierarchy model, can 
be applied for any grand unified model with a type I seesaw mechanism. It ensures that 
the neutrino oscillation data for the mixing matrix will again be satisfactorily fit but now 
with an inverted neutrino mass spectrum. Of crucial importance is the fact that the Dirac 
neutrino mass matrix is the same for both hierarchies, for with a family unification group 
such as SO (10), this mass matrix is closely related to that for the up quark sector, since the 
same set of Higgs representations are involved. This means that the two different hierarchies 
arise solely from the different right-handed Majorana mass matrix textures. Hence the heavy 
right-handed neutrino spectra will differ, as will the possibilities for successful leptogenesis. 




(8) 



Ml = C/;^diag(mi, m'2, m'3)C/t 



(9) 



= -MnM'^'M^. 



Mr = -MIM'-^Mm. 



(10) 
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The nature of the hght neutrino spectrum of course has a direct bearing on the possible 
observation of neutrino-less double beta decay. The relevant parameter is the effective mass 
which can be written as 



(mee) = I (UMNS^fij 



(11) 



j 



The same relation applies for both the normal and inverted hierarchies, with rrij for the 
normal hierarchy replaced by the primed counterparts, m'p in the inverted case, with j — 
1, 2, 3 taken in the same order. 

To make the distinctions more obvious, we shall first use as an example a much studied 
SO{10) grand unified model which explains well the quark mass and mixing data, as well as 
the lepton mass and mixing data with a natural normal neutrino mass hierarchy [8]. Based 
on the 5*0(10) family symmetry with a U{1) x Z2 x Z2 flavor symmetry, one adjoint 45//, 
two pairs of 16//, 16// , along with several Higgs flclds in the 10// and singlet representations, 
the nine quark and charged lepton masses and four CKM quark mixing angles and CP phase 
are well fitted with eight model parameters. The lopsided nature of the down quark and 
charged lepton mass matrices, arising from an electroweak symmetry-breaking Higgs field 
in the 16// representation, readily explains the small Vch quark mixing and near maximal 
[7^3 atmospheric neutrino mixing for any reasonable right-handed Majorana neutrino mass 
matrix, M/j [10]. The type of solar neutrino solution is found to be controlled mainly by 
Mr in this model. In fact, the uncertain nature of this Mr matrix has closely paralleled the 
uncertainty in the neutrino mixing data mentioned in the introduction. When the model was 
initially proposed, a very simple form of Mr led directly to the small mixing angle (SMA) 
solar neutrino solution favored by experiment at that time. A more complicated structure 
with zero subdeterminant for the 2 — 3 sector was later realized to obtain the large mixing 
angle (LMA) solar neutrino solution [11]. The near degeneracy of the two lightest right- 
handed Majorana neutrinos and the resulting possibility of successful resonant leptogenesis 
was recognized in [12]. 

We shall simply begin with the relevant mass matrices and refer the interested reader to 
the literature cited for more specific details of this model. The two Dirac mass matrices for 
the charged leptons and neutrinos are given by 



Here the Dirac matrices arc written with the convention that the Icft-handcd fields label the 
rows and the left-handed conjugate fields label the columns. All nine quark and charged 
lepton masses, plus the three CKM angles and CP phase, are well fitted with the eight input 
parameters. 




(12) 



mu ^113 GeV, 

a = 1.83, 

S = 0.00946, 

(j) = 119.4°, 



m/3 ~ 1 GeV, 

e = 0.147, 

S' = 0.00827, 

?7 = 6 X 10-^ 



(13) 
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defined at the GUT scale to fit the low scale observables after evolution downward from 
Agut- As initially proposed, the two parameters and S'j^ were set equal to zero, but we 
shall allow them to assume also non-zero, but very small, values which can lead to a more 
nearly satisfactory result for leptogenesis [13]. Such values considered are small enough so 
as not to destroy the good fit with experiment obtained initially. 

The most general form for the right-handed Majorana mass matrix considered in [11] 
which gives the large mixing angle (LMA) solar neutrino solution is 



M, 



R 




her] 

2 



e 
— e 



ar] 
— e 
1 



A 



(14) 



as explained by the structure of the Proggatt-Nielsen diagrams [14]. The same parameter, 
e, appears as in Eqs. (12) and (13) under the assumption of a universal coupling of the 45// 
Higgs field to all the matter fields. For the parameters originally chosen in the model, this 
matrix exhibits a normal hierarchy which departs from a simple geometrical form when not 
all three of a, h and c are equal. For the simple case where Sjq — — the seesaw formula 
of Eq. (4) determines the light neutrino mass matrix to be 



MnM^^MJj 
( ^ 



a—b 



\ 



1 



62- 



2 



a—b (ci-^b) 
\ b—a 



b—a 
1 



ttiI/Ar. 



(15) 



The rather extreme hierarchy of Mr is more than balanced by the normal hierarchy of Mjv 
and MJ^ to yield a rather mild normal hierarchy for M^, for selected values of a, b and c. 

We proceed to illustrate the connection between the inverted and normal versions of the 
model with a numerical example. Let us choose the following parameters which lead to a 
normal hierarchy. 



a = 0.25 + 0.15i, 6 = 1.20 + 0.90i, 

c = 0.25 + 0.25i, Ar = 2.90 x 10^^ GeV, 

Sj^ = -0.65 X 10-^ S'j^ = -1.0 X 10-^ 



(16) 



and were recently found to yield nearly successful baryogenesis with a resonant leptogenesis 
mechanism [13]. 

With this choice of parameters hsted, the right-handed Majorana matrix is given by 



Mr = 



A.5i X 10-12 



-(1.058 + 0.794i) x 10"^ 
V (1.5 + 0.9i) X 10-6 

X 2.9 X 10^^ GeV, 



-(1.058 + 0.794i) x 10"^ (1.5 + 0.9i) x 10"^ \ 



0.0216 
-0.147 



-0.147 
1 



(17) 
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from which the seesaw formula of Eq. (4) yields 



/ (7.03 - 5.55i) X 10-^ 
= - -0.0954 + 0.0753i 
V-(8.16-4.79i) X 10-^ 

X 0.0440 eV. 



-0.0954 + 0.0753i 
0.238 - 0.165i 
0.341 - 0.130i 



-(8.16-4.79i) X 10-5 \ 
0.341 - 0.130i 
1 



J 



(18) 



Some results obtained in the model for this normal hierarchy version are 

Mg = 2.96 X 10l^ Ms ~ Ml = 3.06 x 10« GeV, 
ma = 51.0, m2 = 8.8, mi = 2.8 meV, 
(mee) = 0.44 meV, t^b = 2.2 x 10"^°, 



(19) 



for the heavy right-handed and light left-handed neutrino masses, the effective neutrino-less 
double beta decay mass (mee), and the net baryon number to photon number asymmetry 
ratio for the Universe, rjB- The very small departures of ^at and 5^ from zero serve to make 
the lighter two heavy right-handed neutrinos more nearly degenerate and thus enhance the 
resonant leptogencsis effect. We refer the reader to reference [12] for the explicit neutrino 
mixing angle results which are in excellent agreement with the present data. 

To obtain a corresponding model with inverse hierarchy, we select m'2 = 0.10 eV, which 
lies below the maximum of 0.23 eV allowed @ 95% c.l. for nearly degenerate neutrinos by 
the recent WMAP data in combination with the 2dF Galaxy Redshift Survey and Lyman-a 
forest power spectrum [15]. The procedure elaborated above in Eqs. (8) and (9) with the 
help of the unitary transformation matrix. 



I 0.799 + 0.353i 
0.446 + 0.115i 
V -0.157 + 0.020i 



-0.198 + 0.443i 
0.196 - 0.772i 
0.028 + 0.360i 



-0.0134 + 0.0391A 
0.130 -0.370i 
0.002 - 0.918i / 



(20) 



determined from the normal hierarchy case as in Eq. (6), then leads to the light neutrino 
mass matrix. 



M' 



( -0.416 + 0.467i -0.744 + 0.604i 

- -0.744 + 0.604i 0.565 - 0.346i 

V 0.317 -0.109i 0.098 - 0.078i 

X 0.0836 eV, 



0.317 -0.109i\ 
0.098 - 0.078i 

1 / 



(21) 



in the inverted hierarchy version considered, from which the required Mjj is determined to 
be 



/ (1.397 + 0.318i) X 10-1° 
(1.340 - 0.025z) X 10-^ 
V -(1.028 + 0.058i) X 10"^ 

X 1.46 X 10^^ GeV. 



(1.340 - 0.025i) X 10-^ -(1.028 + 0.058i) x 10-^ \ 
0.0204 - 0.0006i -0.142 + 0.003i 

-0.142 + 0.003i 1 / 



(22) 
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Prom Eq. (8) and these matrices we find 



m'2 = 0.1000, m[ = 0.0996, m'^ = 0.0866 eV, 

= 1.49 X W\ = 3.10 X 10l^ M[ = 8.77 x 10^ GeV, (23) 
{m'J = 0.044 eV, rj'^ = 8.5 x lO'^^. 

While Mr and M]j do not appear to be significantly different, the opposite signs of the 
12, 21 and 13, 31 elements for the two matrices play an important role. Moreover, the 
2-3 sector subdcterminant fails to vanish in the latter case. This accounts for the huge 
hierarchy in the right-handed neutrinos in the inverted case, which departs from the near 
degeneracy of the two lighter right-handed neutrinos in the normal case. As a result, the 
nearly successful resonant leptogenesis solution with the normal hierarchy is completely lost. 
The heavy right-handed neutrinos in the inverted case span such a huge hierarchy that the 
lightest is much lighter than the 10^ — 10^° GeV required for successful leptogenesis without 
a resonant enhancement [16]. 

More striking are the differences in the light left-handed neutrino sector, where the 
mass spectrum is indeed inverted but nearly degenerate by the choice of m'2, as illustrated 
in Fig. 1. While the normal hierarchy M^, matrix has two near texture zeros in the 11, 
13 and 31 elements, the elements of are all rather comparable in magnitude. This is 
arises since the determinant of is three orders of magnitude larger than that for M^,, 
i.e., \'m'^'m'2'm[\ ~ 700|m3m2mi|. Such an inverted spectrum is highly unstable against 
very small changes in the matrix elements of M^. It is also very unstable against radiative 
corrections upon evolution downward from the GUT scale, since two of the nearly degenerate 
masses have approximately the same CP parity [17]. Mainly for these reasons, we consider 
the inverted hierarchy solution in this model to be much less satisfactory than the normal 
hierarchy solution. 

We regard this result to be a general feature of grand unified models with a type I seesaw 
mechanism. Unfortunately there are very few models in the literature as predictive as this 
with which we are prepared to test this prediction. However, we can further illustrate the 
claim by considering the whole class of models for which the charged lepton mass matrix is 
diagonal in flavor space. Recall that for the specific model examined in detail above, the 
charged lepton mass matrix is lopsided in that basis. 

With the charged lepton mass matrix. Ml, diagonal in the flavor basis, the unitary 
transformation, f/ii,, is just the identity matrix. Successful models in this class then require 
that — Umns by Eq. (7), where for our purposes we can drop the two Majorana phases 
without loss of generality. We set Umns equal to that in Eq. (2) with the approximations 



C23 — -523 — l/\/2, S12 — y2/7, C12 — y5/7 and S13 = 0. By selecting the normal hierarchy 
neutrino masses found earlier, and inverting Eq. (6) we flnd for the light neutrino mass 
matrix. 



= t/;^diag(0.0028, 0.0088, 0.051)^7^^ 

(24) 



0.451 0.192 -0.192~ 
0.192 2.904 2.196 | x 0.01 eV. 
-0.192 2.196 2.904 
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This matrix clearly has the texture for a reasonably stable normal mass hierarchy. 

To investigate a corresponding inverted mass hierarchy version, we again select m'2 — 0.10 
eV and find with the same Uy^ = Umns the new light neutrino mass matrix, 



M'^ = U* diag(0.0996, 0.10, 0.0866) C/t 



/ 9.975 0.011 -0.011 \ 

0.011 9.325 -0.665 
\ -0.011 -0.665 9.325 / 



(25) 

X 0.01 eV. 



Note that if the small off-diagonal elements of this matrix were neglected, the neutrino mass 
hierarchy would be normal with the lowest two levels degenerate. It is only the small 23 and 
32 off-diagonal corrections which convert the hierarchy into an inverted one. This matrix is 
thus much more sensitive to small changes in the underlying right-handed Majorana matrix 
structure for each model in the class considered. On the other hand, the dominant diagonal 
elements set the scale for the determinant of the matrix and thus the product of the three 
neutrino masses, which again is about 700 times larger than that for M^,. 

These results should be contrasted with those for models involving a softly broken Lg — 
— Lj- flavor symmetry. There the light neutrino mass matrix has large 12, 13, 21 and 31 
matrix elements with small values for the other elements. Hence the results are less sensitive 
to instabilities resulting from small perturbations away from the desired entries, and in such 
models no evolution from a high mass scale is typically involved. 

In summary, we have shown that to every successful normal neutrino mass hierarchy so- 
lution of a grand unifled model corresponds an inverted hierarchy solution with exactly the 
same MNS mixing matrix. The inverted left-handed neutrino spectrum illustrated is nearly 
degenerate and is highly unstable to small changes in the parameters for the right-handed 
Majorana mass matrix as well as to radiative corrections upon evolution from the grand 
unified scale. This suggests that future observation of an inverted hierarchy would tend 
to disfavor grand unified models based on the conventional type I seesaw mechanism. On 
the other hand, grand unified models with type II or type III seesaw mechanisms, where a 
left-handed Majorana mass matrix can arise from direct or induced Higgs triplet contribu- 
tions, or models based on a conserved Lg — — lepton number would then be favored. 
Conversely, observation of a normal hierarchy would eliminate the conserved L^. — — 
number models in favor of the grand unified models. Successful leptogenesis in the latter 
type models is also more favorable. 
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Figure 1: Normal hierarchy and two examples of inverted hierarchy. The 
mass eigenstates labeled by vi, z/2, z^a have the flavor content indicated in 
Eq. (3). The states vi and V2 for the inverted hierarchy are so close in mass 
that these states appear to overlap for the mass scale chosen. The dashed line 
represents the largest nearly degenerate levels allowed by the combination of 
WMAP, 2dF Galaxy Redshift Survey, and Lyman-a data. 
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